Driving phase slips in a superfluid atom circuit with a rotating weak link 
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We have observed well-defined phase slips between quantized persistent current states around a 
toroidal atomic ( 23 Na) Bose-Einstein condensate. These phase slips are induced by a weak link (a 
localized region of reduced superfluid density) rotated slowly around the ring. This is analogous to 
the behavior of a superconducting loop with a weak link in the presence of an external magnetic 
field. When the weak link is rotated more rapidly, well-defined phase slips no longer occur, and 
vortices enter into the bulk of the condensate. A noteworthy feature of this system is the ability to 
dynamically vary the current-phase relation of the weak link, a feature which is difficult to implement 
in superconducting or superfluid helium circuits. 
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Many quantum fluids exhibit superfluid phenomena, 
including dissipationless flow and persistent circulation. 
Furthermore, a quantum fluid interrupted by a weak link 
(e.g., a Josephson junction) can act as a nonlinear inter- 
ferometer, allowing the creation of highly sensitive de- 
tectors such as SQUID magnetometers [1[ and superfluid 
helium gyroscopes [2|-[7j. Ultracold atomic gases offer 
new opportunities for control and measurement [8], and 
it is now becoming possible to realize ultracold atomic 
"circuits" in configurations analogous to those mentioned 
above [9]. Such devices provide opportunities for refin- 
ing our understanding of flow and dissipation in quantum 
fluids, and may also prove useful for inertial sensing. 

Several experiments with ultracold atomic gases in 
simply-connected geometries have observed a critical ve- 
locity for dissipation when a defect is moved through 
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the system 



Other experiments have observed 
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Josephson effects across a thin barrier 
recent experiments with Bose-Einstein condensates in 
multiply-connected geometries have demonstrated per- 
sistent currents [io|, the stochastic decay of persistent 
currents in a series of quantized steps [21(, and phase 
slips across a stationary weak link [9| . 

In this Letter, we demonstrate deterministic phase 
slips between quantized circulation states, caused by ro- 
tating a weak link around a toroidal (i.e. annular) con- 
densate (see Fig. [1]). In the rotating frame, the response 
of the system is analogous to that of a static supercon- 
ducting ring, with a Josephson junction, in an external 
magnetic field. In this analogy, the quantized circulation 
of the superfluid corresponds to magnetic flux, and the 
rotation of the barrier corresponds to an external mag- 
netic field. Rotation of the barrier causes a supercurrent 
/ to flow through it. This flow is dissipationless below 
some critical current I c , leaving the quantized circulation 
state in the ring unchanged. At a high enough rotation 
rate / exceeds 7 C , and excitations can occur. For some 
circuit parameters, these excitations lead to a change in 
circulation state via a discontinuous jump in the phase 
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of the condensate wavefunction, i.e., a phase slip 221-124 [. 



FIG. 1: Ring condensate and the weak link created in it by 
the repulsive potential of a blue-detuned (A = 532 nm) laser 
field, (a) Schematic showing the attractive optical dipole trap 
formed by a red-detuned (A = 1064 nm) horizontal "sheet" 
beam and a vertical "ring" beam, (b) Geometry of the "bar- 
rier" beam. The small double-ended arrow indicates rapid (2 
kHz) radial scanning of the beam. The larger single-ended 
arrow indicates the slow (azimuthal) rotation at up to 3 Hz. 
(c) in situ absorption image of the ring condensate, viewed 
from above. The arc of the barrier trajectory, 6, was centered 
over the azimuth of lowest density, indicated by the dotted 
line. The field of view is 84 /am square, (d) Sequence of in 
situ images showing the effect of the barrier at successive 60° 
intervals around the ring. The barrier height is ~ 60 % of the 
chemical potential. Each image in (d) is the average of three 
absorption images, with a 93 /im square field of view. 



Increasing the rotation rate again until I exceeds I c will 
cause another phase slip. The sharp change of circulation 
(or magnetic flux) in the ring in response to rotation (or 
magnetic field) is the key to the exceptional sensitivity of 
devices such as SQUID magnetometers [1] and superfluid 
helium gyroscopes Q- 

The superfluid state of the condensate is given by 
x/j(r) — \J p(r)e l ^( r \ where p is the superfluid density, 
and (j) is the phase. In an inertial frame <fi is related 
to the superfluid velocity by v s = 2^;V(/>, where h is 
Planck's constant and m is the atomic mass. The cir- 
culation around any closed path, P, must be quantized 
according to the Bohr-Sommerfeld rule, § p v s ■ dl = n^o, 
where n is the integer winding number of the circulation, 
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and = h/m. In a ring geometry without a barrier, 
this gives rise to a series of metastable persistent current 
states, with the ground state having n = 0. 

When two bulk coherent regions are coupled by a tun- 
nel junction, the current is a 2tt periodic function of the 
phase difference 7 across the junction [1]. Even non- 
tunnel coupling can lead to single- valued and 2tt periodic 
current-phase relations under certain circumstances [25| . 
The weak link in our experiments is not a tunnel junc- 
tion; in the direction of flow it is thick compared to the 
condensate healing length, and is more analogous to a 
superconducting constriction (e.g., Dayem bridge) (26| . 
The current-phase relation in such a geometry is in gen- 
eral not single-valued or 2tt periodic, but can become so 
when the peak height of the barrier potential, Uh, is on 



the order of the condensate chemical potential, /in (27- 
I29I ]. As discussed below, our present data are consis- 
tent with the current-phase relation of an ideal Josephson 
junction in series with a small linear inductance. 

We create condensates of ~ 6 x 10 5 23 Na atoms in a 
toroidal optical dipole trap (see Fig. EH). The radius 
of the ring-shaped potential minimum, which coincides 
with the peak atomic density, is determined by absorp- 
tion imaging to be R = 19.2(3) fim [30]. The radial 
Thomas-Fermi half-width of the annulus is w = 10.5(8) 
fim. We estimate the vertical half-width to be « 2 fim 
using the measured vertical trap frequency of 600(5) Hz 
and a calculated chemical potential jao/h ~ 2 kHz. The 
azimuthal variation of the trap depth is < 5% of /in, 
resulting in the nearly uniform density profile shown in 
Fig. Die). The condensate is initially formed in a non- 
rotating state (3l| . 

The rotating weak link was created by a focused, blue- 
detuned laser field, which caused the atoms to experi- 
ence a localized (optical dipole) force repelling them from 
the beam. The focus was a spot ~ 8 /im in diameter 
(FWHM), which was moved along a chosen trajectory 
in the plane of the ring, using a 2-axis acousto-optic de- 
flector. Rapid scanning of the spot allows us to create 
effective time-averaged optical potentials [32|. In this 
experiment we created a wide, flat barrier by scanning 
the radial position of the spot across the condensate as 
a triangle- wave, with a scan amplitude greater than 2w 
(see Fig. [TJd). We rotated this time-averaged barrier az- 
imuthally for 1.5 s at constant Q/2tt < 3 Hz, which 
is < 1/10 the angular frequency of sound propagating 
around the ring. During the first 0.5 s, the strength of 
the barrier potential was increased linearly to a value 
Ub ~ Mo/2. It was then held constant for 0.5 s, and finally 
ramped off over the last 0.5 s. To minimize systematic ef- 
fects due to variations in the effective I c around the ring, 
we determined the region of lowest density and centered 
the barrier trajectory over this region (Fig. H(c)). 

The critical current through the weak link, 7 C , is 
reached when the superfluid velocity becomes equal to 
the critical velocity 0, (HJ . Increasing the barrier height 
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FIG. 2: Quantization of circulation in the ring condensate, 
(a-f) Absorption images of the condensate after 10 ms time- 
of- flight expansion, showing examples of different circulation 
states, n. The field of view is 180 /xm square. The central 
hole is the signature of the persistent current phase singular- 
ity, which is trapped in the center of the ring prior to release. 
(e,f) The off-center holes are due to vortices within the annu- 
lus. (g) Histogram of hole size distribution, showing grouping 
around discrete values due to quantization of circulation in 
the condensate. The vertical axis is the number of occurrences 
per bin, the horizontal axis is the area of the hole [371 ] - This 
histogram shows the same data set used to generate Fig. [3] 



reduces the local superfluid density, decreasing the criti- 
cal current. The ability to dynamically vary the strength 
of the weak link and its current-phase relation is an ad- 
vantageous feature of this atom circuit, potentially allow- 
ing for "third-terminal" functionality that is difficult to 
achieve in superconducting circuits [3J-l36|. 

To determine the circulation state of the condensate, 
we released it from the trap, allowing it to expand in 
time-of- flight (TOF). We then imaged the condensate, 
looking perpendicular to the plane of the ring. A con- 
densate in a non-circulating state expands into a smooth 
density profile that is peaked at the center, as shown in 
Fig.E^a). 

Circulation in the condensate manifests itself in two 
different ways. A pure persistent current state exhibits 
azimuthal flow around the ring, with no vortex core 
present in the annulus. In TOF, this results in a single 
central hole with a size that increases with the wind- 
ing number of the current around the ring, as shown 
in Fig. (2fb-d). At high rotation rates, vortices are ex- 
pected [38|,l39| to enter the annular region and can remain 
there, appearing as off-center holes in the TOF density 
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FIG. 3: Phase slips between quantized circulation states in 
the ring condensate. The vertical axis is the average winding 
number, (n), as a function of barrier angular frequency, 
for two barrier heights: Ub = 0.50 /io (blue squares), and Ub 
= 0.57 fio (red circles, vertically offset by 1 unit for clarity). 
The solid symbols show (n), as determined from the central 
hole size. The open symbols are (n) , plus the average number 
of annular vortices. The number of runs averaged for each 
data point is at least 20; the error bars are an estimate of the 
la confidence interval |43j]. The grey lines show the positions 
at which phase slips are expected for /3l = 2 (lower) and 1 
(vertically offset), with Q /2ir — 1.26 Hz. These values of Pl 
were chosen to match the — >• 1 steps. 



profile. These are distinguishable from the central hole 
caused by a persistent current, as shown in Fig. [2fe,f). 

The histogram in Fig. [2ff) shows the central hole size 
distribution [40| for the data in Fig. [3j comprising 727 
experimental runs. The quantization of circulation is evi- 
dent from the peaking of the distribution about a series of 
values, as has been previously noted (U (HJ. The width 
of the peaks is due to variation in the TOF expansion 
dynamics, primarily caused by shot-to-shot atom num- 
ber fluctuations. Annular vortices are not stable in the 
absence of the rotating barrier, and leave the annulus on 
a time scale that depends on the condensate temperature, 
the trap geometry, and trap smoothness [4l|, . Under 
the conditions of our experiment, the measured lifetime 
of these annular vortices is about 3 seconds, therefore 
an annular vortex formed during our 1.5 second stirring 
procedure is likely to be present during TOF imaging. 

Using the TOF images, we can determine the change 
in circulation caused by the barrier rotating at a given 
angular frequency Q. Figure [3] shows the response of 
the condensate to the rotating barrier for two barrier 
heights. The lower data set is for a Ub = 0.50(7) /io, the 
upper data set is for a barrier height 14(2)% higher. For 
each Ub and we repeated the experiment at least 20 
times, and averaged the measured winding numbers n. 
For small ft the condensate is unperturbed, and n = 
after the stirring procedure. As Vt increases, the average 
winding number (n) changes from to 1. The value of 



FIG. 4: Normalized current I / I c (a) and circulation energy 
E/Ej (b) from equations Q and j4j), as a function of the 
normalized rotation speed Q/Qo (Qo = fr/mR 2 ), for /3l — 2. 
For this /3, Ej/h ~ 30 kHz. The solid black line is the dy- 
namically stable region for n — 0. The dotted lines show the 
unstable regions. The solid blue lines are the stable solutions 
for n — ±1. The red arrows indicate the expected evolution 
of a system initially in the n — state as Q is increased from 
up to a critical point Q c at which a (dissipative) phase slip 
to n = 1 occurs. 



(n) remains 1 up to a higher angular frequency, where a 
transition to (n) = 2 occurs. A comparison of the two 
data sets shows that by varying the barrier height Ub we 
can adjust the critical current I c . In particular, a lower 
Ub corresponds to a higher 7 C , i.e. a higher critical fl. 

To model our experiment, we consider a ID system in 
a frame co-rotating with the weak link. In this rotating 
frame, the superfluid velocity is v£ = v s — O x f, and the 
Bohr-Sommerfeld quantization condition is 



27rn = 2-7T 



•2tt- 



LrI 



•7(J)- 
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Here = 2ttQR 2 is the external "circulation flux" 

due to rotation at angular velocity ^, Lr = l/pm is the 
kinetic/hydrodynamic inductance of the ring, excluding 
the junction 0,135 I is the circumference of the ring mi- 
nus the effective length of the junction, and /?id is the 
mass per unit length around the ring, which is assumed 
constant outside the barrier region. / is the mass cur- 
rent in the rotating frame. The terms on the right of 
Eq. (pp) are, respectively, the Sagnac phase due to the ro- 
tating frame, the change in phase due to current around 
the ring, and the phase drop across the weak link. The 
last term is determined by the current-phase relation of 
the weak link, which we approximate as an ideal sinu- 
soidal current-phase relation plus a small linear kinetic 
inductance Lb 27-29|: 
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We combine equations (pQ) and (j2j) and obtain 



1 

2^ 



J-c 



sin I — 



(3) 



The parameter /3l = Lh/Lj is the ratio of the total 
hydrodynamic inductance L# = L# + to the (fluid) 
Josephson inductance of the weak link Lj = ^o/27r/ c . 
Figure SJa) is a normalized plot of J(O) as implicitly de- 
fined by Eqn. ([3]), with /3l = 2, which is in the hysteretic 
regime [45|, and approximates the conditions of our ex- 
periment for Uh = 0.5 /io- The circulation energy, E, can 
be written as 461: 



E - l B f 1 



1 — COS 27T 
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where = K^l c j2i\ is the Josephson energy of the weak 
link. Equations ((3]) and (jlj implicitly define i£(f2), plot- 
ted in Fig. HJb), with /3l = 2. For each value of the 
winding number n, the stable branches (solid lines) con- 
tinue up to I = I c . Above this point, the system becomes 
dynamically unstable and dissipation can occur (dashed 
lines). At this point in the plots of Fig.HJ there is a state 
of different n with lower total energy, and a phase slip 
occurs, taking the system to that state. The critical an- 
gular frequency at which a phase slip occurs to n ± 1 
is found by setting I/I c = ±1 in Eqn. ((3]), and solving 
for tt: 



9t 
n 



= n± 



2tt 



(5) 



where VLq = ko/ZkR 2 . The solid grey lines in Fig. [3] show 
the values of for n = 0, 1, 2, for /3l = 1 (vertically 
offset) and /3l = 2 (lower), both with Qo/2tt = 1.26 Hz. 
These values of f3 were chosen to match the observed 
phase slips, and are consistent with calculated estimates 
of Lh and Lj. As seen, decreasing increases I c and 
Pl, causing Q c to increase. 

While the data show a clear signature of a phase slip 
from n = to n = 1, the response of the condensate is 
different at higher Q. We attribute this to the fact that 
we have a wide (w « R) annulus. At higher ^, the veloc- 
ity mismatch between the irrotational flow (v oc 1/r) at 
the inner and outer edges and the rotating frame (v oc r) 
makes it energetically favorable for annular vortices to 
appear [38|, |39| . This is consistent with what we observe 
for Q/2tt > 2 Hz, as shown in Fig. [2fe,f). The presence 
of annular vortices allows the winding number to depend 
on the path P. In particular, the difference in wind- 
ing numbers at the inner and outer edges of the annulus 
equals the number of vortices in the annulus. This effect 
is reflected in the difference between the open and closed 
symbols of Fig. [3l The appearance of annular vortices 
is analogous to the penetration of flux lines into a type- 
II superconductor [39|, . The non- negligible width of 



the annulus also necessitates a correction to the effective 
radius of the ring, and therefore to Qq. For a parabolic 
radial density profile with the dimensions of our trap, 
the correction factor is 1.06 [48]. This gives an effective 
^o/2tt = 1.26 Hz. 

In conclusion, we have deterministically driven single 
phase slips between circulation states in an atom circuit 
with a rotating weak link. This circuit is an atomic ana- 
log of a single-junction DC SQUID. The behavior of the 
circuit is consistent with a model assuming a 2tt periodic 
current-phase relation for the weak link. Furthermore, by 
varying the barrier height, we have shown that the junc- 
tion critical current can be controlled experimentally. 

Other phenomena characteristic of superconducting 
circuits should also be observable in atomic circuits. 
Since rotation in our system is analogous to magnetic 
field in a SQUID, our device represents a proof-of princi- 
ple BEC rotation sensor. Additional interesting examples 
include hysteresis, sensitivity to thermal effects, Shapiro 
steps, and multi-junction devices. Our ability to dynam- 
ically change the current-phase relation may even lead to 
new functionality in controlling atom circuits. 
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